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Summary. We have developed a scattering matrix approach to coherent trans- 
port through an adiabatically driven conductor based on photon-assisted processes. 
To describe the energy exchange with the pumping fields we expand the Floquet 
scattering matrix up to linear order in driving frequency. 
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1 From an Internal Response to a Quantum Pump Effect 

The possibility to vary several parameters at the same frequency but with 
different phases [1] of a coherent (mesoscopic) system opens up new prospects 
for the investigation of dynamical quantum transport. The adiabatic variation 
of parameters is of particular interest since at small frequencies the conduc- 
tor stays close to an equilibrium state: the opening of inelastic conduction 
channels is avoided and quantum mechanical phase coherence is preserved to 
the fullest extend possible. The relevant physics has a simple and transparent 
explanation within the scattering matrix approach. 

The variation of parameters leads to a dynamic scattering geometry. Quite 
generally we consider changes in the scattering geometry as an internal re- 
sponse [2] in contrast to the external response generated by voltages applied 
to the contacts of the conductor, see Fig. 1. In general a linear response con- 
sists both of a response to an external potential oscillations, and response to 
internal potentials. The internal response can be expressed with the help of 
the emissivity v{a,r). The emissivity v{a, r) is the portion of the density of 
states at r of carriers that will exit the conductor through contact a. The 
emissivity relates the amplitude I a {uj) of the current in lead a to the am- 
plitude U(r,u>) of a small and slowly oscillating internal potential. At zero 
temperature we have [2] 
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Fig. 1. External and internal response: An ac current with amplitude I(u>), say, 
in the left lead can arise as a response either to an oscillating potential V(t) — 
V{u)e~ l '^ t + V(— a;)e"* ;t at one of the external reservoirs, or as a response to an 
oscillating potential profile &U(r,t) = 8u(r) (U{uj)e~ lu}t + U(— o;)e 1 "') inside the 
mesoscopic sample. 
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Here e is an electron charge and i 2 = — 1. The integral in (1) runs over the 
region in which the potential deviates from its equilibrium value (typically 
the volume occupied by the scatterer). 

The emissivity is expressed in terms of the scattering matrix S (of the 
stationary scatterer) and its functional derivative with respect to the internal 
potential variation W{r,t) = bu{r) ([/(w)e~ itJt + U(-aj)e iiAjt ) : 
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The scattering matrix is evaluated at the Fermi energy E = fi. The summation 
runs over all the leads (for simplicity here assumed to be single channel) 
connecting the sample to external 7V r = 1, 2, 3.. reservoirs. 

Applying the inverse Fourier transformation to (1) gives the current I a (t) 
flowing in response to a time-dependent internal potential U(r, t). Eq. (1) can 
easily be generalized to find the response to an arbitrary field or parametric 
variation of the scattering geometry [3] . To arrive at a general expression we 
remember that the scattering matrix depends on the internal potential U. 
That in turn makes S time-dependent, S(t) = S[f/(t)]. Thus alternatively we 
can express Eq. (1) in the form 
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Originally in [2] the potential U(r, t) is the self-consistent Coulomb po- 
tential. However the form of (3) tells us that the current I a (t) can arise in 
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response to a slow variation of any quantity (parameter) which affects the 
scattering properties of a mesoscopic sample. For instance the current gener- 
ated by a slowly varying vector potential [4] permits to derive the Landauer 
dc-conductance from Eq. (3). But Eq. (3) is not limited to the linear response 
regime. The current I a (t) given by Eq. (3) is a nonlinear functional of the 
scattering matrix. Therefore, the mesoscopic system can exhibit an internal 
rectification effect, i.e., an oscillating internal potential (or any appropriate 
oscillating parameter) can result in a dc current Idc- Since the elements of 
the scattering matrix are quantum mechanical amplitudes the dc-current is 
the result of a quantum rectification process. This is a quantum pump effect 
[3, 4, 5, 6, 7, 8, 9]. An approach to quantum pumping based on (1) and (3) 
was put forth by Brouwer [3]. 

Since the time derivative enters (3), quantum rectification will work only 
under special conditions [3] . Let one or several parameters affecting the scat- 
tering properties of a mesoscopic sample change adiabatically and periodically 
in time. Then the scattering matrix changes periodically as well. Consider the 
point representing the scattering matrix in the space of all scattering matrices. 
During the completion of one time period this point will move along a closed 
line C. Then the dc current Idea, which is the current I a (t) averaged over the 
time period T = 2tt/uj 7 can be represented as a contour integral in the above 
mentioned abstract space [9]: 



lew 

dc ' Q = 4^ J c 



h^) aa - (4) 



The dc current Idc,a, is non-zero if and only if the line C encloses a non- 
vanishing area T . The easy way to see this is to consider a two parameter 
space with parameters being S and S^. Since S and depend on the same set 
of parameters, S = S({Xj}), = S^({Xj}), j = 1,2,..., N p , (which includes, 
for instance, the shape of a sample, the internal potential, the magnetic field, 
the temperature, the pressure, the Fermi energy, etc.) then to get the cycle 
with T 7^ it is necessary to have at least two parameters X\ (t) = X\ cos(ut+ 
(fi) and X2(i) = X 2 cos(a;t + ^2) varying with the same frequency 10 and with 
a phase lag AiyS = ip\ — ip 2 7^ 0. In particular, if the oscillating amplitudes 
are small, i.e., if the scattering matrix changes only a little across T, then the 
pumped current is proportional to the square of the cycle area [3] : 

_ e^sin(Ay)XiX 2 ^ ^ / dS* aj3 dS a/3 \ 

2. ^\d Xl axJ Xi=0iX2=0 ' (5) 

The very simple and compact expression (4) allows to find the pumped current 
for a wide range of situations. Illustrative examples can be found in Ref. [10]. 
In the following we will now discuss the pumping process from the point of 
view of photon-assisted transport through a mesoscopic system. We will show 
how the interlay between photon-assisted transport and quantum mechanical 
interference results in a quantum pump effect [11, 12]. 
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2 Quantum Coherent Pumping: A Simple Picture 

By nature, the quantum pump effect is a rectification effect. A single pa- 
rameter variation only leads to an ac-current. In a two parameter variation 
the modulation of the scatterer due to the second parameter rectifies the ac- 
currents generated by the first parameter. Rectification is only achieved if the 
driven system can scatter electrons in an asymmetric way. Here the asymme- 
try means that the probability T a p for an electron to pass through the sample, 
say, from the lead j3 to the lead a and the probability Tp a to transit the scat- 
terer in the reverse direction differ, T Q/ 3 ^ Tg a . Then the flux of electrons 
entering the scatterer through lead a and the flux of electrons scattered and 
leaving the system through the same lead a differ from each other, resulting in 
a net electron flow in lead a. The resulting current can be viewed as a result of 
an asymmetrical redistribution of incoming flows between the outgoing leads. 

To clarify the physical mechanism which can lead to asymmetric scatter- 
ing we now emphasize the essential difference between the driven scatterer 
and a stationary one. The key difference is the possibility of photon-assisted 
transport. In the stationary case if an electron with energy E enters the phase 
coherent system then it leaves the system with the same energy E. In con- 
trast, in the driven case while traversing the system an electron can absorb 
(or emit) energy quanta nfioj and thus it can leave the system with an energy 
E n = E + nhu). 

It is important that the electron changes its energy interacting with a 
system which is modulated deterministically. As a consequence the inelas- 
tic processes is coherent. If there are several possibilities for transmission 
through the system absorbing or emitting the same energy (say, nhoj) the 
corresponding quantum-mechanical amplitudes will interfere. Such an inter- 
ference of photon-assisted amplitudes can lead to directional asymmetry of 
electron propagation through a driven mesoscopic sample. 

To illustrate this process we consider a simple but generic example. It is a 
system consisting of two regions with oscillating potentials V\ (t) = 2V cos(wt+ 
Lpi) and V2(t) = 21/cos(u;<: + (^2) separated by the distance L. For the sake 
of simplicity we assume that both potentials oscillate with the same small 
amplitude 2V. Consider an electron with energy E incident on the system. In 
leading order in the oscillating amplitudes only absorption/emission of a single 
energy quantum fko needs to be taken into account. So, there are only three 
scenarios to traverse the system. In the first case, an electron does not change 
its energy, the outgoing energy is i?( out ) = E. In the second case, it absorbs 
one energy quantum, 

£(out) = E + hcj.In the third case, it emits an energy 
Hlj, £(° ut ) = E — Tiuj. Since all these processes correspond to different final 
states (which differ in energy £K out ) from each other) then the full probability 
T to pass through the system is a sum of three contributions: 



T = T (0) (E; E) + T (+) {E + hw;E) + T^\E - huj; E). (6) 
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Fig. 2. For the process in which a carrier with energy E traverses two oscillating 
potentials Vi(t) and V^(f), a distance L apart, and absorbs/emits an energy quantum 
hu> there are two interfering alternatives. The modulation quantum can be either 
absorbed/emitted at the first barrier or at the second one. 



Here the first argument is an outgoing electron energy while the second argu- 
ment is an incoming energy. 

The probability T^°\ like the tranmission probability of stationary scat- 
terer, is insensitive to the propagation direction. In contrast T^ + ' and T^' 
are directionally sensitive. Therefore we concentrate on the last two. 

Let us consider T^ + \ There are two possibilities to pass through the sys- 
tem and to absorb an energy Kcu, see, Fig. 2. The first possibility is to absorb 
an energy due to the oscillation of V\{t). The second possibility is to absorb 
an energy interacting with V2(t). In these two processes an electron has the 
same initial state and the same final state. Therefore, we can not distinguish 
between these two possibilities and according to quantum mechanics to calcu- 
late the corresponding probability we first have to add up the corresponding 
amplitudes and only then take the square. Let us denote the amplitude corre- 
sponding to the propagation through the system with absorbing Huj at Vj as 
A^' + ', j = 1, 2. Then the corresponding full probability is: 



T (+) 



2,+) 



2 



(7) 



Each amplitude (either A^' + ' or A^ 2 ' + ^) is a product of two terms, the am- 
plitude A*- llcc) {E) = e ikL (here k = V2mE/Ti is an electron wave number) of 
free propagation in between the potentials and the amplitude A^ describing 
the absorption of an energy quantum fko at the potential Vj. The amplitude 

_4j + ' is proportional to the corresponding Fourier coefficient of Vj(i). The 

proportionality constant is denoted by a. Therefore we have A^ — aVje~ lip > . 

The probability for an electron going from the left to the right is denoted 
by The probability corresponding to the reverse direction - by Tjj~\ 

Our aim is to show that 

t1 +) ^t1 +) - (8) 
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First we consider Tty. Scattering from the left to the right an electron first 
meets the potential Vi(t) and only then the potential V^t). Therefore if an 
electron absorbs the energy Hui at the first potential it traverses the remaining 
part of the system with enhanced energy E + \ = E + Huj. The corresponding 
amplitude is A^'+l = A[ +) A {frcc) (E +1 ). While if an electron absorbed Huj at 
the second potential then it goes through the system with the initial energy E. 
The quantum mechanical amplitude corresponding to such a process reads: 
A {2,+) = ^ iree \E)A { 2 +) . If the energy quantum is much smaller then the 
electron energy, Tiuj <C E, then we can expand the phase factor corresponding 
to free propagation with enhanced energy E + \ to first order in the driving 
frequency: k + iL w (k + ^) L, here v = hk/m is an electron velocity. Thus we 
have: 

A^ « aVe-^e^+v^, 

(9) 

A {2,+) = e ikL aVe -i<p 2 

Substituting these amplitudes into (7) we obtain the probability to pass 
through the system from the left to the right with the absorption of an energy 
quantum Tilo: 

Tl+) =2a 2 T/ 2 |l + cos^ 1 -^ 2 -^J. (10) 

Now we consider the probability Going from the right to the left an 

electron first meets the potential V^(f) and then the potential V\(t). Therefore, 
the corresponding amplitudes are: 

A {i_,+) =c ikL aVe -i Vl ^ 

At' +) « a Ve-' 1 ^e i ( k+ -) L 
Using (7) and (11) we find: 

Tl+)=2 a 2 t/ 2 |l + cos^ 1 -^ 2 + ^|. (12) 

Comparing (10) and (12) we see that the transmission probability depends on 
the direction of electron propagation as we announced in (8). 

Let us characterize the asymmetry in transmission probability by the dif- 
ference Z\T<+) = Ti+) - Tl+): 

AT^ = Aa 2 V 2 sin(A^) sin (^j . (13) 

In our simple case the emission leads to the same asymmetry in the photon- 
assisted transmission probability: AT^ - ) = AT^ + \ Therefore if there are the 
same electron flows Iq coming from the left and from the right, then a net 
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current I = I (AT^ - ) + AT"( + )) is generated. We assume a positive current 
to be directed to the right. 

We see that the induced current / depends separately on Aip = ipi — 
p 2 and on the factor ujL/v. This is an additional dynamical phase due to 
absorption of an energy quantum fiui. Such a separation of phase factors can 
be interpreted in the following way. The presence of the phase lag Aip (by 
modulo 2n) between the oscillating potentials V\ (t) and V2 (t) breaks the time 
reversal invariance of a problem and hence potentially permits the existence 
of a steady particle flow. The second term emphasizes a spatial asymmetry 
of a model consisting of two inequivalent oscillating regions separated by a 
distance L, and tells us that the interference of photon-assisted amplitudes 
is the mechanism inducing electron flow. Neither a phase lag nor a photon- 
assisted process can separately lead to a dc current. 

The simple model presented here reproduces several generic properties of 
a periodically driven quantum system. First, a driven spatially asymmetric 
system can pump a current between external electron reservoirs to which this 
system is coupled. Second, the pumped current is periodic in the phase lag 
between the driving parameters. Third, at small driving frequency, u> — > 0, 
the current generated is linear in w. In contrast the oscillatory dependence on 
oj of (13) is a special property of the simple resonant tunneling structure [13]. 

3 Beyond the Frozen Scatterer Approximation: 
Instantaneous currents 

In experiments the external electrical circuit to which the pump is connected 
is important [14, 15]. Consequently, it is of interest to understand the work- 
ings of a quantum pump connected to contacts which are not at equilibrium 
but support dc-voltages and ac-potcntials. We now present a number of re- 
sults of a rigorous calculation of dynamically generated currents within the 
scattering matrix approach for spinless non-interacting particles [16]. These 
results permit the investigation of pumping also in experimentally more real- 
istic non-ideal situations. 

We generalize the approach of [2] to the case of strong periodic driving 
when many photon processes are of importance. To take them into account 
we use the Floquet scattering matrix whose elements SF, a p(E n , E) are the 
quantum mechanical amplitudes (times ^k n /k, where fc n = \/2mE n /H, E n — 
E+nhw) for scattering of an electron with energy E from lead (3 to lead a with 
absorption (n > 0) or emission (n < 0) of |n| energy quanta Uw. Like in [2] we 
deal with low frequencies (adiabatic driving) and calculate the current linear 
in to. In this limit we can expand the Floquet scattering matrix in powers of 

U). 

To zcro-th order the Floquet sub-matrices Sf (E n ,E) are merely the ma- 
trices of Fourier coefficients S n of the stationary scattering matrix with time- 
dependent (pump)-parameters, S(t) = 5({Xj(t)}). The matrix S(t) is the 
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frozen scattering matrix, in the sense that it describes the time moment t 
and hence stationary scattering. If all the reservoirs are kept at the same 
conditions (potential, temperature..) then the knowledge of only the frozen 
scattering matrix is sufficient to calculate the current flowing through the 
scatterer, see (3). Under more general conditions knowledge of the frozen 
scattering matrix is not sufficient. We stress that the frozen scattering matrix 
does not describe the scattering of electrons by a time-dependent scatterer: 
only the Floquet scattering matrix does. Since wc are interesting in a current 
linear in u> then, in general, it is necessary to know the Floquet scattering 
matrix with the same accuracy. Thus we have to go beyond the frozen scat- 
tering matrix approximation and to take into account the corrections of order 
uj. As we illustrated in Sec. 2 such corrections are due to interference between 
photon-assisted amplitudes. 
We use the following ansatz: 

S F (E n , E) = S a {E) + V^L^El + hLuAn{ E) + O(co 2 ). (14) 

The matrix A(E,t) [with Fourier coefficients A n (E)] introduced here is a key 
ingredient. As we will see this matrix reflects the asymmetry in scattering 
from one lead to the other and back. The unitarity condition for the Floquet 
scattering matrix leads to the following equation for the matrix A(E, t): 

hw(S\E,t)A(E,t)+At(E,t)S(E,t)) =^P{St;S}, (15) 
where V is the Poisson bracket with respect to energy and time 

nst;S} = 1 ^— --— -j. 

The matrix A can not be expressed in terms of the frozen scattering matrix S(t) 
and it has to be calculated (like S itself) in each particular case. Nevertheless 
there are several advantages in using (14). 

First, the matrix A has a much smaller number of elements than the Flo- 
quet scattering matrix. The matrix A depends on only one energy, E, and, 
therefore, it has A r x A r elements like the stationary scattering matrix S. In 
contrast, the Floquet scattering matrix Sf depends on two energies, E and 
E n , and, therefore, it has (2n max + l) xJV r xJV, relevant elements. Here n max is 
the maximum number of energy quanta Tiu absorbed/emitted by an electron 
interacting with the scatterer which we have to take into account to correctly 
describe the scattering process. For small amplitude driving we have n max w 1, 
whereas if the parameters vary with a large amplitude then n max 3> 1. 

Second, the Floquet scattering matrix has no definite symmetry with re- 
spect to a magnetic field H reversal. In contrast both the frozen scattering 
matrix S and A have. The analysis of the micro-reversibility of the equations 
of motion gives the following symmetry: 
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S(-H) = S T (H), 

(16) 

A(-ff) = -A T (H), 

where the upper index "T" denotes the transposition. In the absence of a 
magnetic field, H — 0, the matrix A is antisymmetric in lead indices, A a p = 

Next, using the adiabatic expansion, (14), we calculate the full time- 
dependent current I a {t) flowing in lead a as follows [16]: 



oo ( 

I a (t) = / d££ U If 0,0 - kA \S al3 (E,t) 

I 



2 



(17) 



e 9t 
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dI af3 (E.t) 
dE 



Here /o, Q is the Fermi distribution function for electrons in reservoir a. We 
assume a current in a lead to be positive if it is directed from the scatterer 
to the corresponding reservoir. Equation (17) generalizes (3) to the case with 
external reservoirs being nonidentical (e.g., having different chemical poten- 
tials, temperatures, etc.). The three parts in the curly brackets on the RHS of 
(17) can be interpreted as follows: The first part defines the currents injected 
from the external reservoirs. It depends on the time-dependent conductance 
G a p(t) — ^\S a p(t)\ 2 of the frozen scatterer and hence it describes a classical 
rectification contribution to the dc current Idea- The second part defines the 
current generated by the oscillating charge Q(t) of the scatterer: 

oo 

ow-EE/d%«^. as) 

a 

Here dN a p/dE is the global partial density of states for a frozen scatterer: 

~d~E~ = 47r {-dE- Sa ^ ~ S ^~dE- ) • (19) 

Apparently this part gives no contribution to the dc current. The third part 
describes the currents generated by the oscillating scatterer. The ability to 
generate these ac currents differentiates a non-stationary dynamical scatterer 
from a merely frozen scatterer. 

The instantaneous spectral currents dI a p/AE pushed by the oscillating 
scatterer from lead to lead a read: 



dJ, 



OC0 



dE 



e - (2huRe[S* a/3 A af3 ] + l -V{S ap] S* af} }) . (20) 



The two terms in this equation have different symmetry properties with re- 
spect to the interchange of lead indices. That is most evident in the absence 
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Fig. 3. Interacting with an oscillating potential the electron system gains modu- 
lation quanta of energy. Absorption of an energy nhui leads to the creation of a 
non-equilibrium quasi-electron-hole pair. The quasi-electron (black circle) and hole 
(open circle) can leave the scattering region through different leads. This leads to 
current pulses of different signs in the corresponding leads. 



of a magnetic field, H = 0. In this case the first term on the RHS of (20) 
is antisymmetric in lead indices while the second term has the symmetry of 
the stationary scattering matrix S and is thus symmetric. Therefore, the ma- 
trix A is responsible for the directional asymmetry of dynamically generated 
currents: 

dig/3 , dl 0a 

dE * dE • [Z1) 
We remark, that if one calculates the dc current generated in the par- 
ticular case when all the reservoirs are at same potentials and temperatures 
(fo,a = fo,p, V a, (3) then (17) (after averaging over a pump period) generates 
Brouwer's result, (4). In this case the matrix A plays no role. Different contri- 
butions of this matrix can be combined as in (15). In contrast, the matrix A 
is important in less symmetrical situations, when the electron flows arriving 
at the scatterer from different leads are different. 

The spectral currents dI a p/dE are subject to the following conservation 
law: 



^ dI a0 (E,t) 



^ dE 



0. (22) 



Such a property supports the point of view that these currents arise inside 
the scatterer. They are generated by the non-stationary scatterer without any 
external current source. The appearance of currents subject to the conserva- 
tion law (22) can be easily illustrated within the quasi-particle picture [12]. 
Let all the reservoirs have the same chemical potential [i a = /i, a = 1, . . . , N T . 
We introduce quasi-particles: the quasi-electrons corresponding to filled states 
with energy E > jj, and holes corresponding to empty states with energy 
E < jU. Then at zero temperature there are no incoming quasi-particles. In 
other words, from each lead the vacuum of quasi-particles is falling upon the 
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scatterer. Interacting with the oscillating scatterer, the system of (real) elec- 
trons can gain, say, n energy quanta Tioj. In the quasi-particle picture this 
process corresponds to the creation of a quasi-electron-hole pair with energy 
nTiLo. The pair dissolves and the quasi-particles are scattered separately to the 
same or different leads, see, Fig. 3. If the scattering matrix depends on energy 
then the quasi-electron and hole are scattered, on average, into different leads 
since they have different energies. Suppose the electron leaves the scattering 
region through lead a and the hole leaves through lead (3. Since electrons and 
holes have opposite charge the current pulses created in the leads a and (3 
have different sign. As a result a current pulse arises between the a and /3 
reservoirs. In this picture it is evident that there is no incoming current and 
the sum of outgoing currents does satisfy the conservation law (22). 

Note, from (15) and (20) it is obvious that a conductor with strictly energy 
independent scattering matrix does not produce current and thus it does not 
show a quantum pump effect. 

The current I a {t), calculated with (17), satisfies the continuity equation: 

a 

and thus conserves charge. To demonstrate this we use the unitarity of the 
frozen scattering matrix, J2 a \S a f}\ 2 = l^ap] 2 = 1, and the definition of 
the charge Q(t) of the scatterer, see (18). 

It follows from (22) that the dynamically generated currents dI a p/&E 
do not contribute to (23). Therefore they have nothing to do with charg- 
ing/discharging of a scatterer. The existence of these currents is an intrinsic 
property of dynamical scatterer. 

In conclusion, we clarified the role played by photon-assisted processes 
in adiabatic electron transport through a periodically driven mesoscopic sys- 
tem. The interference between the corresponding photon-assisted amplitudes 
makes the transmission probability dependent on the electron transmission 
direction in striking contrast with stationary scattering. To consider prop- 
erly this effect we introduced an adiabatic (in powers of a driving frequency) 
expansion of the Floquet scattering matrix and demonstrated that already 
linear in u> terms exhibit the required asymmetry. 

The ability to generate currents is only one interesting aspect of quantum 
pumps. Recent works point to the possibility of dynamical controlled gen- 
eration of entangled electron-hole states [17, 18]. This brings into focus the 
dynamic quantum state of pumps. This and other properties will likely assure 
a continuing lively interest in quantum pumping. 
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